Quantum memory for microwave photons in an inhomogeneously broadened spin 
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We propose a multi-mode quantum memory protocol able to store the quantum state of the field 
in a microwave resonator into an ensemble of electronic spins. The stored information is protected 
against inhomogeneous broadening of the spin ensemble by spin-echo techniques resulting in memory 
times orders of magnitude longer than previously achieved. By calculating the evolution of the first 
and second moments of the spin-cavity system variables for realistic experimental parameters, we 
show that a memory based on NV center spins in diamond can store a qubit encoded on the |0) and 
|1) Fock states of the field with 80% fidelity. 

PACS numbers: 03.67.Lx, 42.50.Ct, 42.50.Pq 
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Ensembles of electronic spins have been proposed as 
quantum memories in hybrid architectures for quan- 
tum computing including superconducting qubits \mtM- 
Progress in this direction was reported in a number of 
experiments, demonstrating first strong coupling of an 
ensemble of spins in a crystal to a superconducting res- 
onator 0"11|, and more recently reversible storage of a 



single microwave photon in the spin ensemble 12, [l3j |. 
From these results it clearly appears that inhomogeneous 
broadening of the spin ensemble is a major obstacle, 
which needs to be overcome for hybrid quantum cir- 
cuits to fully benefit from the long spin- coherence times. 
Due to inhomogeneous broadening, quantum information 
leaks from the " bright" collective degree of freedom cou- 
pled to the cavity into dark modes of the spin ensemble 
[14l4l6| . An appealing possibility is to actively and co- 
herently restore it using refocusin g te chniques, inspired 
from magnetic-resonance methods [17| and based on the 
application of ir pulses to the spins acting as time re- 
versal. However, these ideas face a number of challenges: 
(i) The spatial inhomogeneity of the microwave resonator 
field may make it difficult to apply a tt pulse efficiently to 
each spin, (ii) after the 7r-pulse inversion, the spin ensem- 
ble should remain stable despite its coupling to the cavity, 
and (iii) the whole statistics of the collective spin must 
be restored at the single quantum level. The present 
work proposes a protocol, which addresses all these is- 
sues, and we exemplify its feasibility for the specific case 
of NV centers in diamond [HI , using currently available 
experimental techniques and realistic parameters. The 
proposed memory extends the storage times by several 
orders of magnitude compared to 12 , [l3| ■ It is intrin- 
sically multi-mode and thus allows to store reversibly a 
number of quantum states, paving the way to the real- 
ization of a genuine quantum Turing machine 0, [lij . 

In our proposal the ir pulses are performed by rapid 
adiabatic passage [20| through the electron spin reso- 



nance, a method known to tolerate an inhomogeneous 
microwave field. Stability of the ensemble after inver- 
sion is ensured provided the cavity quality factor is suffi- 
ciently low [21| . Since this is incompatible with a faithful 
transfer of quantum information from the cavity into the 
spins, we propose to use a cavity with a quality factor 
that can be tuned in-between the steps of the protocol, 
as was recently demonstrated with SQUIDs 22]. In ad- 
dition, inspired by a recent proposal of atomic-ensemble 
quantum memories for optical photons (23j . we employ 
two 7T pulses in the refocusing scheme. To avoid emitting 
a microwave echo from the inverted spin ensemble, which 
would otherwise be more noisy than the original quantum 
state |24| , we detune the cavity from the spins in-between 



the two pulses (effectively " silencing" this noisy first echo 
[23I] ) . The second echo, formed in a non-inverted ensem- 
ble, restores quantum information with a fidelity up to 
80 % for realistic parameters. 

The proposed physical setup is shown schematically 
in Figs. jUa) and (b); a diamond crystal containing NV 
centers [18j | is placed on top of a transmission- wave-guide 
cavity whose frequency co c [22l | and coupling to the mea- 
suring line k [25l | can be tuned on a nanosecond timescale 
using control lines (not shown in Fig. [1}. The crystal is 
subjected to a constant bias magnetic field Bnv, lift- 
ing the degeneracy of the mg = ±1 states, and bring- 
ing the — > 1 transition to an average frequency uj s — 
2ir ■ 2.9 GHz. In the frame rotating at uj s the free evolu- 
tion of the cavity field and the spin ensemble is described 
by the Hamiltonian: Hq = A cs aja c + Y]j ^-a^z \ where 
a c is the cavity field annihilation operator, A cs = uj c — uj s 
is the (adjustable) spin-cavity detuning, Aj = luj — lu s , 

(1) 

ujj the resonance frequency of the jth spin and a z the 
corresponding Pauli operator. NV centers are coupled 
by hyperfine interaction to the nuclear spin of their ni- 
trogen atom (having a spin 1), causing the mj = — > 1 
transition to split into a triplet separated by Ahf s /27r = 
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FIG. 1. (Color online) (a) Quantum memory circuit. The 
resonator, with frequency cj c and damping rate k tunable at 
the nanosecond scale, is coupled to the spin ensemble (fre- 
quency oj s ) with an ensemble coupling constant g emi . Drive 
pulses of amplitude f3(t) are applied to the spins via the res- 
onator, which can be initialized in a well-defined quantum 
state \a). (b) Amplitude of the microwave field generated by 
a coplanar resonator with quality factor Q = 100 and driven 
by a pulse of 100 [iW power. A static magnetic field Bnv 
is applied parallel to the spins, which are distributed uni- 
formly throughout the crystal, (c) Sub-ensemble distribution 
/(A) of spin-resonance frequencies (circles) consisting of three 
hyperfine-split Lorentzian lines. The solid line shows the ex- 
citation probability for the chosen secant hyperbolic inversion 
pulses (see text), (d) (Solid line) Calculated coupling-strength 
distribution function p(g)g 2 . (Histogram and circles) Sub- 
ensemble distribution used in the calculation. The low- and 
high-frequency cut-offs in p(g) originate from, respectively, 
high (40 (im) and low (0.5 urn) cut-offs in the distance from 
the resonator to NV centers. 



2.2 MHz [26j. In addition, they are coupled by dipolar 
interactions to a bath of magnetic dipol es 1271. which 
is known to govern their coherence time [28|-[30|. This 
bath broadens each of the hyperfine resonances, with a 
Lorentzian line shape [29j | of width w, corresponding to 
a free- induction-decay time T 2 = — . A Hahn-echo pulse 
sequence [13] partially refocuses this coherence, yielding 
a coherence time T 2 which can be several orders of mag- 
nitude longer than T 2 . In this work, we thus model 
the system by the static inhomogeneous spin distribu- 
tion shown in Fig. Hfc) of characteristic width r ~ w 
[27l | . and damped at a rate 71. = T^T 1 in the Markov 



approximation 27]. The spin-cavity interaction is de- 
scribed by: Hi — Y^j <7j(o+ &c + <5"_ where the cou- 
pling constant gj of the jth spin is distributed as shown in 
Fig. Old) [27} . This distribution is of no concern for stor- 
ing the quantum state [l2[; however, it prevents the ap- 
plication of a "hard" it pulse since each spin has a differ- 
ent Rabi frequency for a give n drive amplitude. So-called 
hyperbolic secant pulses [3l|, where the pulse amplitude 
and phase are modulated as a c = a™ ax [sech(/3 soc h£)] 1+v , 



are known to remedy this issue [32|. The pulses are ap- 
plied by an external drive ft modeled by the Hamiltonian 
H cxt = iV2n(l3al — (3*a c ). Note that to achieve the de- 
sired temporal dependence of a c , j3 must be further tai- 
lored in order to account for the cavity filtering and the 
coupling to the spins [27 1. 

The quantum memory protocol, shown schematically 
in Fig. [H[a), aims to store a cavity-field state given at 
t = Q and retrieve it again at t = T mcm with the cavity 
tuned to a "target frequency" A cs . This quantum state 
could be delivered by, e.g., a super-conducting trans- 
ition qubit along the lines of [l2j]. The cavity state is 
then transfered to the spins by setting A cs = for 
a time T swap after which the cavity is parked at Ap.. 
In a lowest-order approximation T swap = ir/2g ens where 
<?cns = [/ g 2 p(g)dg] 1 ^ 2 corres pon ds to the resonator-spin 
ensemble swap rate [13, [l3l Il6|. but in reality is opti- 
mized numerically. For a high-fidelity storage, we set 
« = Kmin = 9^ — with Q max = 10 4 so that the spin 
ensemble and resonator are in strong coupling. Next, in 
order to refocus the reversible spin dephasing we apply 
two 7r pulses at ~ T "^"" and ~ 3T ^ em with A cs = 0, and to 
stabilize the inverted spin ensemble, we set k = K max = 
2 q c . with Qmin = 100 before the ir pulses so that the 

cooperativity parameter fulfills C = < 1 [2l| . An ad- 
ditional constraint comes from the fact that tuning the 
cavity frequency or quality factor with SQUIDs is possi- 
ble only if the cavity field is sufficiently low (|(a c )| ^ 10), 
which requires sufficient delay to allow it to decay after 
the 7T pulses. Between the two ir pulses, we set A cs = A P S 
in order to silence the first spin echo [23j . After the sec- 
ond 7r pulse the quantum state is retrieved from the spin 
ensemble by setting A cs = during T swap after which the 
cavity is tuned to Aji s . 

The numerical calculation of the dynamical evolution 
is made tractable by dividing the spins into M sub- 
ensembles along the lines of [21| keeping account of the 
mean values and covariances between cavity-field quadra- 
tures, X c and P c , and spin components, Sx , Sy m \ and 
Sz of the mth sub-ensemble [27| • Such a representation 
is convenient for determining the memory performance 
for, e.g., coherent input states. Specific for our NV-center 
example we use gens = 27r • 3.5 MHz, w = 2ir ■ 2 MHz cor- 
responding to T* = 0.16 us, T 2 = 100 us [13], and hy- 
perbolic secant ir pulses truncated at a duration of 1 us 
with fi, = 3.5 and /x/3 sec h = 27r • 7.5 MHz. We assume that 
a microwave drive of peak power up to 100 |a.W can be 
applied to the sample input without causing too much 
heating. 

Typical results of our calculations are shown in Fig. [2] 
Panel (b) shows the mean values of X c and P c when a 
weak coherent cavity-field state is given at t = 0. Even 
though the cavity field is very strong during the inver- 
sion pulses at t ~ 2.5 [is and t ~ 7.5 (as, it relaxes to 
negligible levels prior to memory retrieval. Due to an 
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FIG. 2. (Color online) (a) Schematic timing of pulses and 
cavity parameters, A cs and k. Periods of resonance (A cs = 0) 
are marked by gray areas, (b) Cavity-field mean values, X c 
(black) and P c (gray) versus time. The inset re-plots the 
dashed-line region with | (a c ) | on the logarithmic vertical scale, 
(c) The g-weighted transverse-spin-component mean value 
Sf = ^/S^ 2 + Sf 2 (black) normalized to N, the excita- 
tion probability p cxc (gray, dashed curve), and the g- weighted 
excitation probability p^xc (g ra yi solid curve). 



imperfect storage process [marked by the arrow in panel 
(b)] a minor part of the field is left in the cavity (14 
% in field strength or 2 % in energy units), but most 
importantly \{a c )\ recovers at t = T mcm a value compa- 
rable to the one at t = 0. Regarding the spin state, 
we consider the effective, (/-weighted spin observables, 
^ ff = Sif-^ ) V = x,y and g = g C ns/VN, which 
couple directly to the cavity field a c through the interac- 
tion Hamiltonian H\. Panel (c) shows the magnitude of 
these transverse spin components; in the storage part it 
grows as the quantum state is swapped from the cav- 
ity and then decays within T 2 * due to inhomogeneous 
broadening. Despite the excitation of very large mean 
spin components by the it pulses, the much weaker mean 
values of the stored spin states are recovered as a pri- 
mary echo [arrow in panel (c)] and at the final memory 
retrieval. Panel (c) also shows the excitation probabil- 
ity Pexc — S *2N N ano - the effective, g-weighted excitation 

probability pf®. = ^(Ej + N "> versus time - The 

latter reaches 89 % between inversion pulses and levels 
off at 8 % after the second inversion pulse. 

The above results can be extracted from mean-value 
equations alone and demonstrate the feasibility of the 



spin ensemble as a classical memory. In order to assess 
the quantum properties of the memory we also calcu- 
late the evolution of variances by the coupled first- and 
second-moment equations detailed in the Supplementary 
Material 27 [, see Fig. [3] Panel (a) shows the summed 
variance of X c and P c , which deviates from the coherent- 
state value of unity during inversion pulses. At the mem- 
ory retrieval the variance also increases when the cavity is 
tuned to resonance with Q = Q max due to emission from 
spins left in the excited state by a non-perfect inversion 
process (in analogy to Ref. 24]), but most importantly 
this excess noise of only 11 % maintains easily the quan- 
tum nature of the memory. Panel (b) shows the summed 
variance of the spin components, S° and Sy, which re- 
laxes almost to the coherent-state value at the memory 
retrieval. 

We stress the indispensable role of the spin-frequency 
inhomogeneity, which for a resonant cavity in low-Q 
mode gives rise to the effective coopcrativity parameter 

2 i 1 

C = - cna r « 0.38. According to [21( this ensures (i) that 

the excess variance of X c , P c , S% s , and Sy S converge to 
moderate, finite values during the resonant, inverted pe- 
riod (see, e.g., Fig. |3][a,b) at 3 (is < t < 4 \is), and (ii) 
that mean values of the coupled spin-cavity system ob- 
servables relax sufficiently fast from possibly imperfect ir 
pulses as exemplified in the inset of Fig. [U[b) ■ For the 
off-resonant cavity the first and second moments of the 
spin components are damped on the T 2 * time scale as seen 
in Fig. [2jc) prior to the primary echo and in Fig. [3fb) 
at t s» 4 (is, respectively. This is essential for the perfor- 
mance of the protocol; any reminiscence of the inversion 
pulses and excess noise in the spin ensemble must vanish 
both at the time of the primary echo and of the memory 
retrieval. 

To assess the performance of the quantum memory, we 
repeat the above simulation with various other coherent 
input states. A selection of these are shown in Fig. ^c) 
in terms of retrieved mean values and variances (gray 
circles) as compared to those of the input states (black 
circles) . We confirm that the input / output relations con- 
stitute a linear map, which (i) essentially maps vacuum to 
vacuum (with a slightly increased variance) demonstrat- 
ing that the remains of the inversion pulses are negligible 
and (ii) which presents a gain factor Q = 0.79 for the 
mean values. The quadrature variances of the retrieved 
states amount to 2<r 2 = (5X^) + (<5P C 2 ) = 1.11. Since any 
quantum state can be expressed as a superposition of 
coherent states the memory should work for arbitrary in- 
put states, e.g. Schrodinger cats [33|, and qubit states 
encoded in the |0) and |1) Fock states of the cavity. The 
storage time depends on the quantum state and the de- 
sired fidelity. Following [34[ we obtain a qubit fidelity 
F q = 80 % for T mcm = lO iis. 

To investigate the implications of the limited peak 
power available for inversion pulses, the above-mentioned 
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FIG. 3. (Color online) (a) Summed variance of cavity-field 
quadratures, (b) The summed, (/-weighted spin-component 
variance normalized to 2N. (c) Various input states (black) 
and output states (gray, sign reversed) examined in the pro- 
tocol. The center of circles mark mean values whereas the 
radii mark the standard deviation a of the state, (d) Open 
symbols: The dependence of gain Q (blue triangles), qubit 
fidelity _F q (magenta circles), effective excitation probabil- 
ity p°xc ( re d squares), and summed variance 2a 2 (green di- 
amonds) on the peak power P pca k of the external driving field 
during inversion pulses. Closed symbols: Simulations with 
-Ppcak = 100 u.W and homogeneous coupling g — 2% ■ 12.5 Hz, 
leading to Q = 0.82, 2a 2 = 1.02, and F q = 87 %. 



analysis is repeated for a selection of peak powers ranging 
from 20 uW to 500 uW leading to the results presented in 
Fig. |3td) with open symbols. Furthermore, a simulation 
is carried out at P pca k = 100 uW but with a homogeneous 
distribution of coupling strengths, g/2n = 12.5 Hz (solid 
symbols in Fig. GUd)). Clearly, increasing Pp ea k presents 
an increase in Q due to a better inversion process, but 
since in an intermediate regime a fraction of spins expe- 
riences a poor inversion process due to insufficient Rabi 
frequency (limiting the inversion performance illustrated 
by the dashed curve in Fig. [2jc)) we observe the non- 
monotonous behavior of 2<r 2 shown in Fig. ^d). While 
increasing driving powers may be infeasible from an ex- 
perimental point of view an alternative route to improve- 
ment lies in tailoring a more homogeneous distribution of 
coupling strengths, e.g. by limiting the distance between 
NV centers and the cavity. 

Continuing the analysis with a homogeneous coupling- 
strength distribution (solid symbols in Fig. E^d), F q = 
87%), we find the limiting factors for the obtained fidelity, 
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FIG. 4. (Color online) The cavity field \a c \ versus time in 
a multi-mode storage example with four input fields (|o c | = 
3, 0, 1, 2) separated by 0.29 (is, with memory time 12 u.s. The 
amplitude-cross-talk is below 3 %. 



which in terms of gain can be written approximately as: 
g = Go expt-K^ + 2T chirp ]) exp(-7 ± [T mom - 0.7 \is\). 
The re-dependent factor yields s» 0.92 due to cavity de- 
cay during the resonant swapping process and during the 
initial and final frequency chirp of duration T c hi rp - The 
7±-dependent factor yields ~ 0.91 due to spin decoher- 
ence (partly suppressed when the quantum state resides 
in the cavity or a population degree of freedom). The 
main contribution to excess noise arises from imperfect 
inversion processes, e.g. due to the dephasing rate 7^ 
during ir pulses. In the limit T2,Qmax 00 the qubit 
fidelity becomes w 97 %, and the origin of the remaining 
infidelity [Go « 0.97 and 2cr 2 rj 1.01) includes a non- 
perfect cavity-to-spin transfer (arrow in Fig. HJb)) and 
residual imperfections in the inversion processes. 

As demonstrated experimentally for classical pulses 
[35] ]. the spin-ensemble quantum memory is multi-mode 
in nature, which we confirm by simulating storage and 
retrieval of four pulses, see Fig. |H The number of stor- 
age modes (proportional to T% /T 2 ) that can be faithfully 
addressed and refocused is estimated to be ~ 100 [27j |. 

In summary, a multi-mode spin-ensemble-based quan- 
tum memory for cavity fields has been proposed and an- 
alyzed for a specific realization using NV centers in di- 
amond. With realistic experimental parameters a qubit 
fidelity of F q = 80 % is predicted for T mom = 10 us. The 
main limiting processes are clarified, and we predict a 



qubit fidelity Fq > |, better than achieved by any clas- 
sical strategy, for memory times of T mcul < 69 u.s. The 
memory time may be further increased by dynamical de- 
coupling techniques 36 , 3?| or quantum-state transfer to 
nuclear-magnetic degrees of freedom [38] . 
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The effective spin-i Hamiltonian 

The interaction Hamiltonian between the spin state of NV centers and the microwave 
cavity field is calculated in the following. The electronic part of the Hamiltonian for a single 
NV-center in an external magnetic field reads: 



O- H = DS 2 z +E(S 2 x ~S 2 y )+g NV ^ B S-B, (1) 

(N 



where D = 2n ■ 2.88 GHz is the zero-field splitting between the ms = and ms = ±1 states, 
E is a splitting between the m s = ±1 states induced by non-axial strain g NV = 2 is the 
gyro-magnetic ratio, /iB is the Bohr magneton, «S is the real spin (not to be confused with 
the collective ensemble-spin S defined elsewhere in the manuscript), and B is the external 
magnetic field. The NV-center axis defines the direction of quantization along z, which we 

a '. 

^ ■ assume in this work to be parallel to the sample. 

The NV center is effectively reduced to a two-level system by applying a static bias 
magnetic field B z . If the Zeeman energy verifies g^rv^B-Bz ^ E, which is the case for 
B > 0.5 mT, the spin states = 0, ±1 are also energy eigen states. Taking as our ground 
state |g) = \m s = 0) and excited state |e) = \m s = 1), the transition energy becomes: 
to = D + gNV^BB z + O (^ 9n ^ bB j ■ As will be detailed below, the quantized cavity field is 
linearly polarized within the xy-plane and can be written: B = <5B(a c + aJ). In perturbation 
theory the interaction part of the Hamiltonian then reads: 



Hi — ^nv/^b 

S'NV/^B 



SJB x {r) + S y 8B y {v) 



[a c + at) 



} v S + {5B x {r) - i5B y {v)} + ^{^(r) + iSB y (r)}\ (a c + at c ) 

= 3^ [a c a + {5B x {v) - i5B y {v)} + a\a.{8B x {x) + i5B y {v)}] <2) 
V2 

5f NV /i B |5B ± (r)| , , A , 1 
->■ -j= [a c o + + <a_J . 

The second line uses the standard definition of raising- and lowering operators, S± = <S x ±iS y , 



with properties: S±\ms) = \/<S(<S + 1) — mj(ms ± l)\ms ± 1). Restricted to our effective 
two-level system, |g) and |e), these operators can be stated in terms of Pauli operators: 
S± = v^2o"±, which has been exploited in the third line together with the rotating- wave 
approximation. The last expression is obtained by writing 5B x (r) + i5B y (r) = \5~B±(r)\e ie ^ 

1 



and re-defining our basis states as |g) and e~ ld ^\e) (such that the phase factor is absorbed 
into er±). We conclude that the coupling between the cavity field and the effective spin-| of 
the jth NV center located at rj is given by: 

^NV/iB^B^r,-)! 



9j 



V2 



(3) 



Analytical expressions for the magnetic field of a wave traveling in the z direction along 
a coplanar transmission line of the geometry, shown in Fig. 1(b) of the main text, can be 
expressed as a function of the voltage V between the center conductor and the ground 



plane 
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In this expression, fj, is the vacuum permeability, i] = 376.7 the vacuum impedance, c the 
speed of light, e e fr = (e r + l)/2, e r being the relative dielectric constant of the substrate, 
and S = W/b, 8 = (S + W)/b, F n = ^ are geometrical quantities with S the width of the 
coplanar waveguide central conductor, W its distance to each ground plane [see Fig. 1(b) of 
the main text], and b the ground plane width. 

In a A/2 coplanar waveguide resonator of length L defined by a transmission line open at 
z = and z = L, the fields are generated by voltages and currents V(z) = V(0) cos(7r,z/L) 
and I(z) = —iV(0)sm(7rz/L)/Z , so that the current in the middle of the transmission 
line verifies I(L/2) = iV(L)/Z , the same relation as between the current and voltage of 
a traveling microwave at a given position of a coplanar waveguide. The zero-point rms 
fluctuations of the magnetic field in the middle of the resonator are thus directly obtained 
by replacing in the previous expressions Vqjby the zero-point rms fluctuations of the voltage 
at the resonator end 5Vq(L) = u c \JTiZq/t^ [3]. In this work we chose the waveguide geometry 
so that Zq = 50f2, as was the case in a number of recent experiments 4j. 

Using the expression for the rms zero-point magnetic field fluctuations, we obtain the 
dependence of the coupling constant of a spin located in the middle of the resonator (z = 
L/2) as a function of x, y. From that we calculate the coupling constant distribution shown 
in Fig. 1(d) of the main text. 



The electron-spin-resonance transition of the NV ensemble 



We now discuss the spin lineshape and coherence properties that were chosen in the main 
text. The electronic spin of all NV centers is also coupled by hyperfme interaction to the 
nuclear spin of the nitrogen atom. This leads to an additional structure of the electron spin 
resonance. In the case of NV centers with a 14 N nucleus as is the case with diamond samples 
prepared with natural nitrogen abundance, the spectrum is known to consist of three peaks 
separated by A hfs /27r = 2.2 MHz 

In addition to the hyperfme coupling to the nitrogen nucleus, NV centers are coupled to 
a bath of other spins that determines their coherence properties, namely the line width and 
the spin echo time T^. Two different well-identified spin baths contribute to NV center de- 
coherence: nitrogen impurities (so-called PI centers), and 13 C nuclear spins. Hybrid circuit 
experiments usually require a rather high concentration in NV centers (typically 1 ppm or 
more), in order to efficiently absorb single photons or equivalently to reach the strong cou- 
pling regime. This high concentration also implies a comparatively large concentration of PI 
centers (on the order of a few ppm) , which at this level are mostly responsible for the sample 
coherence properties A sample with [NV] ~ 2 ppm and [PI] ~ 2 ppm can be achieved 
with usual crystal preparation techniques (electron irradiation of a nitrogen-doped crys- 
tal). According to [6], such a sample would have a line width of each of the hyperfine peaks 
around w /2ix = 0.5 MHz, a spin-echo coherence time T2 ~ 100 |os, and according to Q] would 
be coupled to a coplanar resonator with an ensemble coupling constant g ens /2-K ~ 3 MHz. 
These are parameters similar to the ones we use in our calculations (see main text), apart 
for the line width which we take somewhat broader w/2ir = 2 MHz. Broadening the line 
width purposely could easily be achieved, either by adding a slightly inhomogeneous mag- 
netic field, or by a slight misalignment of the magnetic field with a crystalline axis, which 
results in a slight difference in Zeeman shift for NV centers with different orientation [7] . We 
finally note that in order to avoid complications of Electron Spin Echo Envelope Modulation 
(ESEEM) caused by the interaction with the bath of 13 C nuclei [8j, it would be preferable 
to use a 12 C-isotope-enriched diamond crystal for this experiment. With a sample having 
natural abundance of 13 C, our refocusing protocol would still work, but only for discrete 
times corresponding to the spin-echo revivals depending on the applied magnetic field, and 
for NV centers having their axis well aligned with the magnetic field js, 9]. All effects taken 



together, the spin distribution can be expressed as 



/(A) 



w 
Qtt 



(A - A hfs ) 2 + 



+ A 2 + ni + (A + A hf s) 2 + 
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The corresponding characteristic width 



10| can be derived: 
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2 J (7± + f) 2 + |A 2 fs ' 



(6) 



We note that the free induction decay of an ensemble of identically prepared spins, 
(<r_ (0)) = (<t_(0)), can be expressed as: 



N 



(5_(t)> = J3<^ ) (0)>e- 



(7±+iAj)t 



N(&-(0)) / f(A)e- { ^ +iA)t dA 



(7) 



-(5_(0))[l + 2cos(A hfs t)]e-^ + f) t . 



As stated in the main text, the time scale for free induction decay, T 2 * = ^ (taking 7j_ <C w), 
is seen to characterize the envelope function of the modulated decay, caused by the hyperfme 
splitting. 



EQUATIONS OF MOTION FOR FIRST AND SECOND MOMENTS USING SUB- 
ENSEMBLES 

With the considerations of the above section, a Hamiltonian for the physical system can 
be written in the frame rotating at the central spin frequency, u s as: 

n N N 

H = HA cs ala c + - ^ A^ + ihV2K(pa< c - (3*a c ) + 9j{^+a c + <3"i j) aJ), (8) 

3=1 3=1 

where A cs = u c — u s is the detuning of the cavity- resonance frequency co c from u s , and 
Aj = Uj — u s with Uj being the individual spin-resonance frequencies. The sums add the 
contribution from the N spins residing in the volume considered. An external coherent-state 
field, (3, may be used to drive the cavity field a c through a coupling capacitor, which gives 
rise to the cavity-field-decay rate k. The c- number is normalized such that |/3| 2 is the 
incoming number of photons per second. 

Decay mechanisms are treated in the Markovian approximation. In our previous work [lo| 
we employed the formalism of quantum Langevin equations, which has pedagogical appeal 

4 



for understanding the structure of Eq. (fl5|) below. However, the present work describes 
more general states of the spin ensemble, in which case we prefer the approach of the 
Lindblad master equation. This reads: ^ = -^[H,p] + ^ k T^[ck\P with V[c k ]p = — \c\c k p — 
+ c k pc k , where Ci = \p2ka c accounts for the cavity leakage. In addition, for the jth 
spin a population decay with rate 7y is modeled by C2j = ^/7j[<3"^' ) and phase decay with 
characteristic time r by c 3j - = -j^b^ ■ Even though 7y = in the main text, it is included 
here for generality. 

Our aim is to describe both the storage of weak cavity fields in the spin ensemble and 
the impact of very strong driving fields needed for spin-refocusing schemes. To this end 
we develop a formalism, which is applicable for any saturation level of the spins and which 
accounts for the first and second moments of relevant physical operators [Eq. (J2j)]. The 
inhomogeneous spin ensemble is divided into M sub-ensembles, JAi, M.% y . . . , M.Mi which 
can each be regarded as homogeneous with coupling strength g m , spin resonance frequency 
A m , and containing N m spins for m = 1, . . . , M. 

Next, the dynamical variables are described in terms of the real- valued operators: 



a c + a' c 



V2 



p = -i[a c - a' c ) 



V2 ' 



Mr 



(9) 



Mr 



M, 



The X c and P c operators describe the quadratures of the cavity field with [X c , P c ] = i, 
while the components correspond to twice the total spin in each sub-ensemble with 
[«S'- m ^ , S^} = 2itjkiS\ m \ These operators are now linearized around their mean values: 
X c = X c + 5X C , P c = P c + 5P C , and = + for k — x,y, z, i.e. by definition, 

(SX C ) = 0, etc. Using the master equation, the following mean-value equations can then be 
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deduced: 



dX c 
dt 

dPc 
dt 

dt 



dS, 



(m) 



dt 

dSi m) 



dt 



-kX c + A CS P C - + 2v^/3 R , 

-kP c - A CS X C - J2 + 2v^A, 

m v 

- 1± S^ - A m S^ - V2g m (S^P c + (SS^SP C )), 

- 7 ±^ m) + A rn S^ - V2g m (Si m) X c + (8S^8X C )), 

V2g m (Si m) P c + (SS^5P C ) + S^X C + (5S^5X C )) 
- 7|| (^M+iV w ), 



(10) 
(11) 

(12) 
(13) 
(14) 



where /3 R = Re{/3} and A = Im{/3} and j± = ± = ± + We note that Eqs. (fT2> 
( |T4l) contain covariances between the spin and cavity-field operators. In practice, for all 
calculations in the main text, these are very small and can be omitted such that the mean- 
value equations by themselves form a closed set. 

In total, there are 2 + 3M operators to consider (two quadratures of the cavity field 
and three spin components for each sub-ensemble), which can be written on vector form, 
y, with X c and P c as the first two entries and S^, Sy 1 ^, Sz , Sy 2 \ Sz , etc., on the 
following entries. The covariance between any two of these operators, y k and yi, is defined 
as C(y k ,yi) = 2Re{{Sy k Syi)} , and in particular C(y k ,y k ) = 2Vai(y k ). These are grouped 
into a (2 + 3M) x (2 + 3M) covariance matrix 7 for which the entry of the kth row and Zth 
column is C(y k ,yi). The time-evolution of 7 is: 



^ = M 7 + 7 M T + N, 

dt 



(15) 



where the driving matrix M is given by: 
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-v / 2^ m ^ m) 
^f^9mSy ^ y/2g m Si 

-7_l -A m -V2g m P c 
A m -7 ± —y/2g m X c 
V2g m P c y/2g m X c 



The noise matrix N is block diagonal: 



N 



V 

u« 










xj( m ) 



V 



2k 
2k 



27|| Si' 



(m) 



4 7± iV m 

47_LiV m A -, !rv 
2 7 ||^ m) 2 7 ,|4 m) 4 7! |(^ m) +iV„ 



2 7 || 4 m) 



(17) 



(18) 



The above derivation of Eq. ( 1151) follows from the master equation with the approximation 
that moments of order higher than two have been neglected. Since second moments have 
proved negligible in Eqs. f lT2|) - f|T4|) . it can be expected in succession that moments of order 
higher than two are also negligible and the approximation is well justified. 

Note that both M and N contain mean values, and hence the second-moment equation 
do not form a closed set. However, if all spins are close to the ground or excited state, 
~ ±N m , the S^-operators are effectively removed from the dynamics and the first- 
and second-moment equations reduce to the linear, de-coupled sets published in lol |. 



IMPLEMENTING THE QUANTUM MEMORY PROTOCOL 



In the following the timing details of the quantum memory protocol are stated. The 
schematic setup shown in Fig. 2(a) of the main text is presented in further detail in Fig. |5j 
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FIG. 5. The detailed pulse sequence (not to scale) showing the evolution of the external control 
parameters, A cs , k, and (3. T^ v marks the time, which the quantum state spends effectively in the 
cavity. 



The protocol sequence consists of 21 parts, and the duration of the jih part is denoted 
by Tj. There are three important cavity-frequency detunings in the protocol. (1) At the 
"target frequency" A* s = 2rr ■ 100 MHz the cavity field is given to us for storage or we must 
deliver the state at this frequency after memory retrieval, (2) at resonance A cs = the spin- 
cavity interaction is effective, and (3) at AJ? S = 2rr ■ 50 MHz the cavity is parked when the 
spin-cavity interaction needs to be turned off (in safe distance from the target frequency). 
Maintaining a chirp rate of 2rr • 10 MHz/ns, we denote by T£ = 5 ns the time of the chirp 
— > AP g and by = 10 ns the time of the chirp — > A* s (equal to T chirp in the main text). 
Hence, in the protocol we take T x = T 2 \ = and T 3 = T 6 = T 9 = T 13 = T 16 = T 19 = T£. 
The cavity-field decay rate k is also varied linearly in time with a duration taken to be 
^5 = Tiq = T12 = Tn = T K = 10 ns. 

The duration T 2 = T 2 o = T swap = 73.7 ns of the storage and retrieval is determined 
numerically by the condition that the cavity field amplitude \(a c )\ is minimized after part 
3 of the sequence. The "ideal" time of = 71.4 ns is modified due to cavity decay, 
inhomogeneous broadening, and the fact that the information transfer extends slightly into 
parts 1 and 3 (or parts 19 and 21). 

The duration of the external driving pulses are chosen to be T-j = Tn = = 1 \is. The 
desired intra-cavity-field mean value a c is obtained by tailoring the external driving field 
(3 as follows: Solve first Eqs. (}l2"]) - (fli|) with X c and P c given by the desired a c (the secant 
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hyperbolic pulse stated in the main text) and with initial conditions Sx = = and 
g(m) _ where plus and minus are chosen for the second and first it pulse, respectively. 
Second, the external /3 can be found from re-arranging Eqs. ffTUl) and ffTTT) [we write it on 
complex form]: 



*2k 



— + (k + iA cs )a c + % g m Sl 



(19) 



(m) S {m> —iS (m> 8 

where S_ = — — g - 2 — * s the solution from the first step. The first term ^ accounts 



is the solution from the first step. The first t< 
for the finite response time of the cavity while the third term including describes the 
reaction field from the magnetic dipole moments of the spin ensemble. 

After each inversion pulse the mean values of the cavity field quadratures, X c and P c , and 
of the transverse spin components, S x and S y , will inevitably contain small, non-zero values 
due to a non-perfect inversion process. These mean values must relax to appropriate levels 
before the cavity parameters can be changed (due to saturation effects in the cavity-tuning 
circuitry), which is ensured by (parts 8 and 15) during the additional time on resonance, 
T 8 = T 15 = T res = 1 (xs. Moreover, these mean values together with the possible excited-state 
excess variance must be reduced even further: If the stored cavity field was the vacuum state, 
the spin-mode coupled to the cavity must also be effectively in the vacuum state at the time 
of the primary spin echo in part 11 and at the time of the memory retrieval. Otherwise, any 
stored quantum state will be distorted. Specifically, for the second inversion pulse the time 
T15 + . . . + Tig must be long enough for the mean values to decay and the time Ti 6 + . . . + Tig 
must be long enough for the excess noise to decay (a similar condition applies after the first 
inversion pulse). 

Since relaxation is required after both the first and second inversion pulses, it is convenient 
to arrange the sequence such that the primary spin echo occurs right in the middle of the 
inversion pulses. This is obtained when the total time before the first inversion pulse equals 
the duration after the second, T\ + . . . + T 6 — Ti 5 + . . . + T 2 i, which with the choices above 
reduces to: 

T 4 = T rcs + T 18 . (20) 

Next, the duration T echo , between the effective focus points of the reversible spin dephasing 
process, must fulfill (see bottom part of Fig. [5]): 
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Techo T mem 2T cav 2 ^ ^ Tj, (21) 

j=7 



where the last expression is equal to twice the time between the inversion pulses. Finally, 
the total memory time T mem is given by us: 



21 



^mem — ^ Tj . (22) 
3=1 

The above expressions ( 120|) - (!22|) present four equations with the four unknowns T 4 , T n , T 18 , 
and T echo . After some algebra we find: 

T u = - 2TI - 2T K - T n - T res - T°g, (23) 

'J 1 r T 1 T^efE 

rp _ -'mem rpt OTP T 1 n T 1 T 1 i cav /0/A 
J 18 — — T -'a A k 2 Jrcs_ J swap H — , {^) 

and T 4 and T echo follow immediately from Eqs. f )20|) and f[2Tj) . The effective time T^ v can 
be calculated by choosing some initial guess, e.g. halfway into part 2, and extending Ti 8 
such that the time £ rev ivai of the spin revival occurs within part 18. Then the correct value 

follows: T^ v = trcvival ~ ^echo- 

The multi-mode version of the protocol discussed around Fig. 4 in the main text has 
been simulated as follows: After the initial storage (parts 1-3) and a waiting time of 0.19 [is, 
the cavity is chirped to the target frequency A* s and the cavity field is then replaced by 
a new coherent field ((SX^) = (<5P C 2 ) — > |, correlations between the cavity field and each 
sub-ensemble of spins is set to zero, while spin-spin correlations are maintained). This new 
quantum state is then stored in complete analogy to the first input state. The retrieval is just 
the reverse process, and when the cavity reaches the target frequency A* s , the cavity is reset 
to the vacuum state. This protocol models a perfect, infinitely fast swap between the cavity 
and some "target system" , which holds the unknown quantum state in the storage part of 
the protocol and which is prepared in its ground state at the point of retrieval. In addition, 
we used a cavity Q-parameter of 10 3 between these pulses to facilitate a lower cross talk. 
In the simulation shown in Fig. 4 of the main text, the gain factor is 0.80, being lower by 
the factor e~ 7 - L< - 2 ^ due to the extra memory time as compared to the case of homogeneous 
coupling strengths in Fig. 3, and the variance is maintained at 2a 2 = 1.02. As mentioned in 
the summary of the main text, the experimentally feasible case with inhomogeneous coupling 
strength and P pea k = 100 |J.W presents non-classical storage capabilities with T mem up to 69 
(xs, which can be deduced using the fact that an additional memory time AT will decrease 
the gain by the factor e~ AT//T2 . Using the pulse separation of 0.29 \is as in Fig. 4 the protocol 
will accommodate 100 pulses when T mem = 69 [is. 
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